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A Global Jet/Circulation Model for Young Stars 
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Abstract. Powerful, highly collimated jets, surrounded by bipolar molecular outflows, are commonly observed 
near Young Stellar Objects (YSOs). In the usual theoretical picture of star formation, a jet is ejected from a 
magnetized accretion disk, with a molecular outflow being driven either by the jet or by a wider wind coming 
from the disk. Here, we propose an alternative global model for the flows surrounding YSOs. In addition to a 
central accretion-ejection engine driving the jet, the molecular outflow is powered by the infalling matter and 
follows a circulation pattern around the central object without necessarily being entrained by a jet. It is shown 
that the model produces a heated pressure-driven outflow with magneto-centrifugal acceleration and collimation. 
We report solutions for the three different parts of this self-similar model, i.e. the jet, the infalling envelope and 
the circulating matter that eventually forms the molecular outflow. This new picture of the accretion/outflow 
phase provides a possible explanation for several observed properties of YSO outflows. The most relevant ones are 
the presence of high mass molecular outflows around massive protostars, and a realistic fraction (typically 0.1) of 
the accretion flow that goes into the jet. 
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1. Introduction 



Star formation occurs in molecular clouds. Observations 
have shown that the accretion phase, during which the 
central object builds up its mass, is very often accompa- 
nied by the powerful ejection of prominent bipolar out- 



the controversy are the following questions. Has the cen- 
tral jet or wind sufficient thrust to drive the molecular 
outflow? If yes, how is the momentum transferred to the 
CO lobes? 



flows, that are comprised of fast jets ( 


Ray et al. 1996, 


Eisloffcl et al. 2000)|, 


Reipurth et al. 2000 


) and slower and 


less collimated molecular outflows ( 


Bachillcr 1996). A 



There seems to be 
a magneto-centrifugal 
thcr launched from the 



a consensus 
origin, and 
accretion 



Wardle fc Konigl 1993 



topic of great theoretical and observational interest is the 
dynamical connection between the jets and the molecu- 
lar outflows. Strong arguments, such as the large masses 
observed in the CO flows, suggest that most of the molec- 
ular outflow material consists of swept up ambient gas 
(Bachiller & Tafalla 1999). The only possible source for 
this gas is the surrounding ambient cloud. This implies 
that a driving agent accelerates molecular cloud mate- 
rial and gives rise to the observed CO lobes. An ob- 
vious candidate for such an agent is the central fast 
jet. The dynamical relationship between the optical jet 
and the surround ing bipolar molecular outflows, however , 
is controversial (Bncll et al. 1980, Rodriguez et al. 1980, 



(Blandford fc Payne 1982 



Paatz fc Camenzind 1996] , 



that jets have 
they are ei- 
disk ( disk wind) 



Konigl fc Pudritz(2000) 





Pelleticr & Pudritz 1992 




Sauty & Tsinganos 1994 






Ferreira 1997 


Contopoulos & Sauty 2001 



the interaction of the pro- 
with the disk (X-wind) 



Bally fc Lada 1983, Welch et al. 1983). At the heart of 



or from the location of 
tostar's magnetosphere 

( |5hu et al. 1994] , |Ferreira et al. 1999[ , |Shu et al. 200fj| ). In 

the disk wind model, the engine consists of a Keplerian 
disk threaded by a magnetic field that is either generated 
in situ, or advected in from larger scales, while X-winds are 
magnetised stellar winds where the interaction of a proto- 
star's magnetosphere with the surrounding disk results in 
the opening of some magneto-spheric field lines. On the 
other hand, the precise mechanism generating molecular 
outflows is poorly understood. The latter are generally 
believed to be driven by a jet (e.g. Padman et al. 1997, 
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Chernin & Masson 1991) or by a wide-angle wind 



(Barral fe Canto 1981, Smith 1986, 5hu et al. 1991) 



but all the proposed possibilities present difficulties 
(Cabrit & Bcrtout 1992). For example, in the case of a 
wide-angle wind with a steady wind-blown cavity, the 
shocked wind converging at the tip of the cavity should 
form a jet which is not observed. However, in the case of 
the wide-angle X-wind model, Shang, Shu & Glassgold 
(1998), for example, have shown that a Herbig-Haro jet 
appears along the rotation axis, i.e., , in the central 
part of the cavity. On the other hand, for such models, 
the ambient medium has a density that decreases like 
r~ 2 away from the core, which cancels the decrease in 
wind ram pressure due to the radial expansion, ensuring 
a constant shell velocity. This, however, does not seem 
to be consistent with the sha pe observed for the largest 
flows flCabrit fc Bertout 1992; ). 



In the case of jet-driven molecular outflow models, or 
"prompt entrainment" mechanism, the bulk of the molec- 
ular outflow is accelerated ambient gas near the head 
of the jet or more precisely along the wings of its asso- 
ciated bow shock {e.g. Chernin & Masson 1995). This 
idea is supported by some recent simulations (Micono 
et al. 2000), where it is shown that a large fraction 
of the initial jet momentum can be transferred to the 
ambient medium. However, in some cases, the jet may 
fail to provide the necessary momentum, as in L1551 
IRS 5 ( Fridlund fc Liseau 1998| ). Moreover, observations 
frequently suggest that jets decelerate near their outer 
edges, probably due to entrainment of the surrounding 
material, but show little or no deceleration near the 



flow axis (Giovanardi et al. 2000). Also, in some cases, 
the axis of the cavity differs from the axis of the jet, 
e.g. by about 20° for HL Tau (|Close et al. 1997b|), or by 



about 40° in the case of NGC 2261 flClose et al. 1997b 



Canto et al. 1981 ). Another major problem appears to be 
that st eady jets produce molecul ar lobes with large aspect 
ratios ( |Chernin fe Masson 1991 ). This would mean that 
the jet should be changing its position over time to accel- 
erate different parcels of molecular gas ("wandering jet"), 
or varying its flow velocity resulting in internal shocks 
which push material sideways. Finally, based on numer- 
ical simulations, the efficiency of prompt entrainment in 
accelerating molecules is found to be poor in the case of 
YSO jets with densities comparable to their ambient en- 
vironment (Downes & Ray 1999). 

In addition to these problems, many bipolar out- 
flows from massive protostars transport masses largely ex- 
ceeding those of their associated stars. The bipolar out- 
flow masses appear to range from about 10 to 4800 so- 
lar masses, with a mean value around 130 (Churchwell 
1997). This probably indicates that much of the outflow- 
ing molecular gas could not have been accelerated close 
to the central YSO, since the protostar only represents a 
small fraction of the outflowing mass. Most of the mod- 
els, discussed previously, face difficulties in explaining such 
massive outflows. Indeed, in this case, it seems likely that 
deflection of infalling matter into bipolar outflows may 



be a crucial mechanism, since entrainment and swept-up 
mass do not appear to be able to account for the very large 
observed outflow masses. This is the mechanism that we 
propose in the present article. 

Finally, we point out that a correlation exists be- 
tween the bolometric luminosity of the central exciting 
source and the rate of momentum injection into the out- 
flow flBally fc Lada 1983| , pabrit et al. 1997| ). These quan- 
tities are roughly consistent with a single power law 
(with slope of about 0.7-0.8) across the full range of 
source luminosities extending over five orders of magni- 
tude ( pabrit et al. 1997 ). This continuous trend suggests 
that similar mechanisms may be responsible for the pro- 
duction of molecular outflows from both low and high 
mass stars. As we will show, the present model can ac- 
commodate this observational constraint. 

We start by presenting the global model in section || 
and by writing down our set of self-similar MHD equations 
in section Then, in section [|, we show the properties of 
the solutions along the streamlines, while, in section^, we 
discuss the main quantities that characterize the flow such 
as the fluxes of mass, momentum and energy, the plasma 
beta or the electric current. In section ^| we present some 
physical consequences of the model and we discuss sev- 
eral observational implications. Finally, we give our con- 
clusions in section 0. 



2. The Model 

We propose a global model in which the central jet is 
driven by a magneto-centrifugal accretion-ejection engine, 
while the molecular outflow is powered by the infalling 
matter through a quadrupolar circulation pattern around 
the central object. The situation is schematically repre- 
sented in Fig. [|. The molecular outflow may still be af- 
fected by entrainment from the fast jet, but this would 
be limited to the polar regions and it would not be the 
dominant factor for its acceleration. 



circulation model (II) ( 


Hcnriksen & Valls-Gabaud 1994 


Fiege & Henriksen 1996, 


|Lery et al. 1999 


Aburihan, Fiege, Henriksen & Lery 200l|) and an in- 



falling envelope model (III). We suggest that the gas in 
all the three regions I, II, and III can be described by the 
self-similar heated, quadrupolar and axisymmetric magne- 



tohydrodynamic model ( 


Henriksen & Valls-Gabaud 1994 


Fiege & Henriksen 1996, Lery et al. 1999 


), described 



rate the jet model inside the self-similar circulation 
model, and to include the infalling envelope model 
near the equator. The torus and the jet were not 
part of the previously p ublished models of this type 



( Fiege fc Henriksen 1996| , Lery et al. 1999 ), and their 



inclusion is a new result reported here. Therefore a 
description of the flows around YSOs can be obtained 
from the axis to the equator with a single set of equations 
written in the framework of a self-similar hypothesis. 
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Fig. 1. The different parts of the model. Schematic 
representation of streamlines of the flows around YSOs 
in the poloidal plane. The zones respectively correspond 
to (I) the jet, (II) the circulation region, and (III) the 
infalling envelope. The thin accretion disk, the axis as well 
as the innermost region (r < r a « 100 AU) of the protostar 
containing the accretion-ejection engine are represented in 
grey and are not treated in the present model. 

Note that the model does not directly treat the thin 
accretion disk near the equator nor the accretion-ejection 
engine that provides the jet in proximity of the central 
source, since the physical processes governing these re- 
gions are substantially different from the large-scale in- 
fall/outflow region. 

2.1. The Circulation Model 

The c irculation model (II in Fig. [j]) has already been stud- 
ied in Lery et al. 1999, and we only recall its main features 
here. In the infall phase, the flow is slowed down by the 
increasing radial pressure gradients due to heating by the 
central source, and by an increase with proximity to the 
central object of the strength of the barrier due to the 
rotation and the magnetic field. This pressure 'barrier' 
deflects and accelerates, by means of its poloidal gradi- 
ent, much of the infalling matter into an axial outflow. 
The Poynting flux included in the model increases both 
the velocity and collimation of the outflows by helping 
to transport mass and energy from the equatorial to the 
axial regions (see Lery et al. 199£ for details). The out- 
flow velocity is naturally more pronounced as the axis of 
rotation is approached since the pressure gradients are 
strongest there and the density is smallest. These gradi- 
ents also act to collimate the flow, forming a "throat" in 
the axial region. 



The main attraction for the circulation model is that 
large quantities of neutral material can be accelerated di- 
rectly into the outflow without the need for jet entrain- 
ment, whose plausibility has long been debated. Moreover, 
the model can generally explain the observed opening 
angles and velocity structure seen in high mass systems 
( [Ficge fc Hcnriksen 199(| |Lery et al. 19990 . 

Self-similar models cannot be globally consistent in 
time even if they are non-stationary, since they are, in 
general, ignorant of the initial conditions. The self-similar 
circulation model is intermediate in space between a re- 
gion around the symmetry axis containing the jet, and the 
infalling envelope close to the disk. These regions are de- 
scribed in the following. Note that the axis itself and the 
equatorial plane are strictly excluded from the domain of 
self-similarity although they may be approached asymp- 
totically. 



2.2. The Jet Model 

In most jet models, the matter is thought to be transferred 
from an accretion disk and/or a thick torus surrounding 
the central forming star, to the jets. In the most widely 
accepted scenario the outflow acceleration is mediated by 
magnetic fields that thread Keplerian type disks. In the 
present model, no detail is given about the origin of the 
jet itself but our description is compatible with most mod- 
els mentioned in the introduction. We consider that the 
system consisting of the protostar and the thin accretion 
disk acts as a reservoir for the jet and only study the 
flow inside the jet (I in Fig. |l|). The same set of equations 
used for region II can be applied to the jet. As for the 
circulation solutions, the velocity is larger closer to the 
axis but the dynamics is dominated by the variations of 
the various quantities in the radial direction, i.e., along 
the jet. The jet model applies only to a limited range of 
finite radii, where the self-similarity applies. Specifically, 
our model cannot address the launching mechanism of the 
jet near the star. The structure of a supersonic jet far from 
the source is disconnected from its origin. Thus, it is not 
essential to include a detailed description of the central 
engine in our model. 

2.3. The Self-similar Infalling Envelope Model 

During the formation of the star, it is expected that an 
accreting disk forms and feeds the central protostar from 
infalling ambient cloud core material; these disks may have 
radii from 20 AU up to several hundred AU, and masses 
containing a fraction of the central protostar mass. Hence, 
the last element of the model is the region where material 
is falling from the cloud onto the central object. It consists 
of a standard flattened dense accretion disk, not treated 
by the present model (i.e., the equatorial-most region in 
Fig. [j]), and of an extended and less dense torus of infalling 
gas and dust (III in Fig. [l]). The latter can be as large as 
a few thousand AU and is also described by the set of self- 
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similar equations of the circulation model. The material, 
for both the circulation and the infalling envelope models, 
originate from the same region in the molecular cloud. We 
thus suggest that the same set of equations may apply in 
both cases. This statement would not be applicable to the 
thin disk where resistivity, ionization fraction and other 
parameters can substantially differ from the outer regions. 



3. The Equations of the Model 

3.1. The Variables 

The solutions are developed within the context of r-self- 
similarity wherein a power of r multiplies an unknown 
function of 9; the spherical coordinates r, 9 and <j) being 
used. The only physical scales that enter into our calcu- 
lation are the gravitational cons tant G, the fixed central 
mass M, and a fi ducial radius r Q ( Fiege fc Henriksen 1996 , 
Lery et al. 1999|) . Note that Fiege et al. expressed r a in 
terms of G, M, and the central luminosity L*, while Lery 
et al. offered a different choice. In the present work, we 
leave scale that can be freely specified for any 

observed outflow source. We also introduce two other pa- 
rameters that are the limiting angles between the jet, the 
infall and the circulation zones. The first angle, 6 m im goes 
from the axis of rotation to the boundary of the jet (angle 
between zone I and II in Fig. [j]). The second angle, 9 max , 
corresponds to the limit of the circulation region (angle 
between zone II and III in Fig. @). The power laws of the 
self-similar system are determined, up to a single param- 
eter a, if we assume that the local gravitational field is 
dominated by a fixed central mass. In terms of the fidu- 
cial radial distance, r , the self-similar symmetry is sought 
as a function of two scale invariants, r/r and 9, in a 
separated power-law form. The self-similar index a is a 
free parameter of the solution, but must lie in the range 
— 1/2 < a < 1/4, for simultaneous infall and outflow to 
occur. The radiation field can help us to understand the 
physical meaning of a. Indeed, in this context, the self- 
similar behaviour that we impose on the radiation field 
takes the form F r , 



{GM/r f' 2 M/rl(r/r ) af - 2 f(9) 



( Fiege fc Henriksen 1996| , Lery et al. 1999] ). In this equa- 
tion, the index ay is a measure of the radiation energy 
balance as a function of radial distance. When a/ is neg- 
ative, it corresponds to a loss in energy due to the net 
radiation leaving the region, while, for a positive aj, the 
system has an energy input due to the absorption of radia- 
tion. If one supposes that the opacity is predominantly due 
to dust, the index a/ is related to a by aj ~ —2(1/4 — a). 
Consequently, the case a/ = 0, where there is no net flux 
of radiation through the domain of our solution, corre- 
sponds to a = —1/4. For larger (respectively smaller) val- 
ues of a, there is a net lo ss (gain) of energy by radiation 
( Fiege fc Henriksen 1996 ). 

Hence, if we assume that the gravitational potential 
is dominated by the central mass, i.e., self-gravitation is 



negligible, the equations of radiative MHD admit the fol- 
lowing radial scaling relations for the variables 



GM 



-1/2 



u{0), 



2Q-1/2 




y<t>A e ) ' 



GM 2 



2Q-3/2 



kT 



GM 



m^rriH 



0(0), 



rad 



3/2 



M 



(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



where v, p, B, p, T and F ra( i respectively correspond to 
the velocity, the density, the magnetic field, the pressure, 
the temperature and the radiative flux. In these equa- 
tions the microscopic constants are represented by k for 



Boltzmann's constant, 



for the mean atomic weight, 



run for the mass of the hydrogen atom. In the last equa- 
tion, the index af is a measure of the loss (if negative) 
or gain (if positive) in energy by radiation as a func- 
tion of radial distance. The self-similar variable directly 
related to magnetic field y can be divided into poloidal 
and toroidal components. In the present model the two 
components are not equal. They are respectively defined 
by Vp,<t>{6) = Map^cj, j yj Anfi(9) . Consequently the system 
also deals with two different components of the Alfvenic 
Mach number M ap and M a ^ defined by M% J6) = 

v p </>/ (Bp 4,/^ n pj ■ We will refer to 9 as the value of the 
self-similar temperature Q{0) on the axis. 



3.2. The Equations 

In order to make the system tractable, we assume axisym- 
metric flow so that d/d<j> = and all flow variables are 
functions only of r and 8. We further restrict ourselves to 
steady models (i.e., d/dt = 0). Magnetic field and stream- 
lines are required to be quadrupolar in the poloidal plane 
for the circulation model (see |Lery et al. 1999 for the spec- 
ification of the boundary conditions). Under these assump- 
tions the self-similar equations are: 



1. Mass flux conservation 
Id 

(1 + 2 a) /i u r + 



sin0 d6 

2. Magnetic flux conservation 



(fjtug sin (9) = 0, 



5/4) u r 



1 d ( ug sin 9 



sin 9 d9 



Up 



0, 



(7) 



(8) 
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3. Radial component of momentum equation 



ue 



du r 

w 



M~ p 2 ) - (uj 



1/4 



M 2 



4- 1-— 



0, (9 



4. ^-component of momentum equation 



1/2 



+ 



M 2 
due 



ue- 



dd 

2 



+ ui cot 



Utj, du^ 
M 2 p d0 



M~f - 1 

a<p 

u 2 dy p 
y p M 2 p d9 



u 



dy<t, 



de 
~dJ 



5. Angular momentum conservation 



Qdfi 
~jld9 



1 du v 
dd 



1 



-cot 1 



M ap M at p 
1 



-1/2 
1 



0. 



1/4 



10 



M ap M ai p 



MapMatf, J y^MapMacf, dd 

6. Faraday's Law plus zero comoving electric field 



d(u<pue) 
dl 



u^Ur+u^ue— In 
at) 



1 



1 

y<t> 



0,(11 



0.(12) 



This set of equations can produce either circulating flows, 
pure infall, or pure outflow. Note that the boundary condi- 
tions guarantee mass conservation in the circulating region 
(II), while in the two other cases, mass fluxes are not con- 
served. Indeed, the source acts as a reservoir for the pure 
outflow (region I) while the accretion disk would act as 
sink for the pure infall (region III). Therefore, any discus- 
sion on net mass infall/outflow can only apply to region I 
and III. 

4. Example of Solution 

From a numerical point of view, we integrate the set of 
equations from a given solution close to the axis up to a 
given angle Q m in ■ This solution near 8 — has finite values 
for V r , B r and the density. The other components of the 
velocity and magnetic field are vanishing, otherwise they 
would show jump discontinuities at the boundaries. The 
last point from the jet model solution can be used as initial 
condition for the integration of the same set of equations 
of the circulation model. However, we have restricted our- 
selves to solutions in which the boundary conditions for 
the circulation region are zero for B^, Bg, V^, Vq, a good 
behaviour for the circulation solution. We note that un- 
der these hypotheses no material crosses the boundaries 
defined by the cones delimiting the three regions, and that 
the solutions are in pressure balance. 

Once the circulation solution is obtained, we integrate 
then from this solution down to a given angle {9 max ) in 
the equatorial region, excluding a small zone correspond- 
ing to the thin accretion disk. We also ensure that B^, 




e 



e 



e 



Fig. 2. YSO Jet solution as a function of the an- 
gle from the axis. Velocities and magnetic field are re- 
spectively given in km s _1 and Gauss. Parameters are 
a = —0.12, 0o = 0.35. Here, in order to dimensional- 
ize the quantities, the mass of the central protostar is set 
to one solar mass, and r is kept fixed at a distance of 10 4 
AU from the source. 



Bg, Vtj,, Ve vanish at the interface between the circulation 
solution and the infalling envelope. The main difference 
between the two solutions is the change of signs of the 
^-components of the velocity and the magnetic field. This 
corresponds to a solution that goes towards the axis in the 
case of the circulation solution, and towards the equator 
for the infalling envelope solution. 

The three components of the velocity and of the mag- 
netic field for a jet solution are shown in FigJ| as a func- 
tion of the angle between the axis and a limiting angle 
(Omin = 0.1 in this case), and at a distance of 10 4 AU, 
the mass of the central protostar being set to one solar 
mass. Radial velocities of several hundred of km s _1 are 
obtained near the axis. There, the radial component of the 
magnetic field also dominates the other component. The 
rapid increase in velocity towards the axis is easily ex- 
plained as follows. From Bernoulli's law, a decrease in den- 
sity along a streamline, and hence in pressure, corresponds 
to an increase in velocity. This is what we observe towards 
the axis, where material encounters a low pressure region. 
Bernoulli's law strictly applies only in the absence of a 
Poynting flux, in which energy is directly transferred be- 
tween the magnetic field and the kinetic energy. However, 
it is clear from our solutions that the Poynting flux is not 
dominant in this region. Note that the flatness of b r in 
the innermost part of the jet is not a general behavior of 
the self-similar solutions but corresponds to the present 
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Fig. 4. 2D iso-contours of the density and velocity 
field. The protostar is located in the center. The density 



ranges from 10 to 10 



the lowest densities corre- 



e 



e 



e 



Fig. 3. Circulation and infalling envelope solutions. 

Solid lines and dashed lines respectively correspond to the 
circulation and the infall solutions. Velocities and mag- 
netic field are respectively given in km s~ x and Gauss. 
The distance from the source is taken to be 10 4 AU. 
Parameters for the solution are a = —0.12, O = 0.35, 
9min — 0.1 and max = 1.55, the transition between the 
solutions being at 9 = 1.2. 



set of parameters. It is noteworthy that the strength of 
the radial component of the field could help to stabilize 
the jet w.r.t. internal instabilities (Lery et al. 2000). In 
most of the jet, the angular momentum and the magnetic 
pressure are negligible with respect to the kinetic energy, 
whereas, close to the interface with region II, the hoop 
stress becomes stronger, due to the increase in B$. The 
change of slope in the variables at 6 ~ 0.07 corresponds to 
the location of the steep gradient in B r and of the maxi- 
mum amplitude of in the jet. If this angle is taken as 
a measure of the broadening of the jet during the propa- 
gation at large scales, we may derive a full opening angle 
of the flow of about 8°, which is in the range of values 



sponding to the darkest regions. The size of the box is 
given in AU. 

neither a general result, nor a boundary condition that we 
impose. Solutions with negative values of the radial veloc- 
ity can also be obtained. The rotation is maximum at the 
turning point but it is still dynamically important closer to 
the disk. The ^-component of the velocity shows that the 
gas is moving towards the axial region in the circulation 
model while it is directed towards the equator for the in- 
falling envelope. Finally, the maximum of Vg corresponds 
to a maximum in the toroidal component of the magnetic 
field B^ which produces a radial "hoop-stress" in the out- 



derived from the observations (Eisloffcl et al. 2000)). 

The solutions in the circulation and the infalling en- 
velope torus region are reported in Fig.||. Following the 
flow, the radial velocity is negative in the infalling region, 
dtp > 7r/4, it changes signs at the turning point (the clos- 
est point to the source) located at 9 tp ~ n/4 in the present 
case. Finally it increases and passes from a few km s _1 to 
more than 10 km s _1 in the axial region. Note that the ra- 
dial velocity vanishes in the equatorial region where it con- 
nects to the thin accretion disk, for this solution. This is 



flow that helps the collimation (Lery & Frank 200C). 

We also present the density iso-contours, with super- 
imposed poloidal velocity vectors, in Fig.|] for the entire 
model in the poloidal plane. The density in the jet is lower 
than the density in the molecular outflow by at least a 
decade in the present example, but this is not a general 
result of the model. In the molecular outflow the den- 
sity increases from the outflow to the infalling region, de- 
creases at the interface and reaches again a maximum near 
the equator. Other examples of density profiles as well as 
variations of the main quantities along a streamline can 
be found in Lery et al. (1999) for the circulation model. 

We finally present the streamlines obtained with the 
solution given above in Fig. ||. The streamlines are pro- 
jected in the poloidal plane and are shown with arrows 
that indicate the direction of the flow. The outermost re- 
gion is dominated by the circulation solution except close 
to the axis and in the equatorial region. On the other hand, 
the accretion-ejection is dominated by the dynamics near 
the protostar. 
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Fig. 5. Streamlines around a one solar mass pro- 
tostar. Streamlines projected in the poloidal plane and 
given by the solution for a — —0.12, 9o = 0.35. The solid 
lines correspond to the circulation and the infall models, 
and the dashed lines to the jet model. The arrows indicate 
the direction of the gas on the streamlines. The typical 
sizes of the solution is 10 4 AU. 



5. General properties of the Global Model 

5.1. The Fluxes 

The empirical relationships between the mass, momen- 
tum, and energy fluxes provide important observational 
constraints on physical models of bipolar outflows. Indeed, 
in our model, the total mass flux M, momentum flux P 
and kinetic energy flux E are respectively given by 



Moc(r/r ) 1+2Q ,Poc(r/r ) 



l/2+2a 



, E oc (r/r a ) 



2a 



(13) 



It follows that models with a < —1/4 and a > —1/4, 
i.e., with loss in energy by radiation through the bound- 
aries, have mass and momentum fluxes dominated by the 
inner and outer regions respectively. Moreover, the energy 
flux is dominated by the interior region if a < (and by 
the exterior for a > 0). 

Let us now consider in more detail the mass fluxes 
in this solution. The total mass flux coming from the 
molecular cloud is given by Mi n = Mi n f a u + M acc , where 
Mi n f a u corresponds to the infalling portion of the cir- 
culation pattern, and M acc to the mass flux in the en- 
velope which is subsequently channeled toward the star 
through the accretion disk. The total mass loss rate reads 
M out = Mmo +Mjet, Mmo and Mjet being the mass loss 
rate of the molecular outflow and of the jet, respectively 
Since the infalling gas is conserved along the path in the 



circulation model, Mi n f a u = Mmo- Thus, the global mass 
flux conservation reads 



M 



jet 



M„ 



M* = / M ac 



(14) 



where M* is the accretion rate onto the star, and / is the 
fraction of the accretion flow that goes into the jet. 

The mass fluxes are calculated as M = 
Jg° ut p{0)V r {9)2-KR 2 d9, where 9 in and 9 out are the 
poloidal angles delimiting each region of the system, and 
R a fiducial distance in the radial direction. We thus 
consider fluxes through sections of a spherical surface 
centered on the star, the distance of which from the 
source is arbitrary as long as the material keep flowing 
in conical volumes (far enough from the disk). By using 
this method, the values that we obtain can be tunable 
within a few orders of magnitude of G 1 / 2 M 3 / 2 r^ 3 ^ 2 . For 
the solution presented here, the mass fluxes in the jet, 
the molecular outflow and the accreting torus are found 
to be, respectively: 



M 



jet 



6 x lQ- 7 M Q yr- 



M infall = M MO « 8.5 x lO^Moyr" 1 , 



M acc w 6 x 10- 6 M Q yr-\ 



(15) 
(16) 
(17) 



The value of the mass loss rate in the jet is of the 
same order of magnitude of the ones found in several 



observational studies, 



e.g. 



M 



(Bacciotti et al. 1999), M, 



3.7 10~ 7 for HH34 
2.4 10~ 7 for DG Tau 
( Bacciotti et al. 2001 ) . Our value for Mmo is also similar 
to the ones estimated obscrvationally for typical molecular 



i jet 



flows dBachillcr 1996| ). We also get / = M jet /M acc « 0.1, 



for this particular model, again in good agreement with 
observations and predicted by other MHD models specifi- 
cally treating the launch of the jet. For example, / r* 0.03 
for self-similar disk-wind models (e.g. Ferreira 1997), and 
/ w 0.3 for X-wind models (e.g. Shu et al. 1994). 

Our model only applies to the YSO environment from 
r Q w 10 2 - 10 3 AU to 10 4 - 10 5 AU, which is outside the 
region where the jet is launched. However, we speculate 
that a quadrupolar jet could form on a smaller scale. In 
order to understand the type of acceleration in the jet, 
it can be useful to look at v = A-Krjv^o/ B a , where rj = 
pv r /B r is the mass to magnetic flux ratio. In the example 
presented here, we have v = 0.03 <C 1. This shows that the 
flow is not accelerated magnetically but centrifugally or 
thermally. This could be do ne through scries of giant flare s 
close to the central object ( Feigelson fc Montmcrlc 1999 ), 
but this is beyond the scope of the present work. 

Finally we point out that in our model the poloidal mo- 
mentum flux, calculated as P = f Q out p(9)V^(0)2-nR 2 d9, 
is about the same in the jet and in the molecular 
outflow. This is als o in agreement with observations 
( Bacciotti et al. 1999 ). In our model, however, there is 
no strong dynamical relationship between the jet and the 
molecular outflow. Thus it seems that the observational 
similarity between the momentum rates in such outflows 
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may have no relationship with the mechanism acceler 
ing the molecular lobes. In a future paper, we plan 
constrain the global model with observations of spec 
sources for which the opening angle as well as the m 
fluxes can be estimated. 

5.2. The Energetics of the Flow along Streamlines 

The acceleration of a magneto-hydrodynamic wind is u 
ally described in terms of the centrifugal force acting ale 
the field lines. Alternatively, the development of the fi 
can be described in terms of energies. At the closest po 
to the source, part of the energy is in the form of magnc 
energy, which gets converted (in part) into kinetic ener 
as we will see. In this section and in the next ones we "\ 
describe and study the variations of the energies, of 
plasma beta and of the electric current along streamlines 
in order to illustrate the physical aspects of the model. We 
will concentrate, in particular, on the circulation model 
that is the corner stone of the global model. 

In analogy to Bernoulli's theorem, it is useful for our 
discussion to define a function Be that corresponds to the 
sum of the specific kinetic energy, the specific potential en- 
ergy (or gravitational energy), the specific enthalpy of the 
gas, and the specific magnetic energy, all divided by the 
local specific gravitational energy to make Be dimension- 
less. In the present case, the gravitational energy per unit 
mass is approximately the square of the local free-fall ve- 
locity. One should note that Be is not a streamline integral 
due to the presence of the magnetic field. Nevertheless, it 
remains useful for discussing the flow of energy between 
the various components. This function of 9 is given by 



Eg 

.2 



E 



H 



E B )/(GM/r) 



Be{9) = (E k 

= (ul +uj + ui)/2 - 1 + 6 + u7(8tt/4/ 2 ), (18) 



where u, and fi respectively represent the dimensionless 
velocity, temperature and density. The variables E^, Eg, 
Eb respectively correspond to the kinetic, gravitational, 
and magnetic energies, while Eh is the enthalpy. The vari- 
able y is related to the magnetic field by B oc u/y. The 
parameter a is the self-similar exponent, M is the mass 
of the central object and r Q is a fiducial radius (See pa- 
per I for details). Note that the full energy equation that 
has Be as a first integral has been given previously (See 
equation 15 in paper I or Henriksen 1997). In the flow, the 
positivity of Be suggests that the gas should reach infinity 
with a net positive energy, although this is not strict since 
the energy is not conserved. On the other hand, when Be 
is negative, the gas should not spontaneously escape to 
infinity. Note that the positivity of Be does not imply a 
lack of conservation of energy, but the excess Be repre- 
sents an energy transfer from hotter regions at smaller 
radii to motion, i.e.. mainly from the enthalpy to the ki- 
netic energy. Of course when proper boundary conditions 
are applied and the full global problem is solved, the total 
energy would be conserved. The net gain in Be from the 
infall region to the outflow region in its various compo- 
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Fig. 6. Energetics of the Flow along a Streamline 

The Bernoulli-like function Be (heavy solid line) is rep- 
resented with its various components, the specific kinetic 
energy Ek (dashed line) , the specific magnetic energy Eb 
(point-dashed line) and the specific enthalpy Ep (dotted 
line) in the circulation model. The peaks take place at the 
turning point. The curvilinear abscissa s along a stream- 
line is given in AU (a = -0.12,6 = 0.5, r a « 100AU, 
M « lM©) 



nents (See Fig. ||) represents the energy transferred to the 
gas when it is closest to the star. 

In Fig. ||, the solutions are for a = — 0.12,O = 0.5, 
r Q w 100AU, M w 1M Q . Note that Be is smaller than 
the kinetic energy in this figure because the gravitational 
energy, not represented here, is negative. Far from the 
central object, the dominating energy is gravitational at 
the beginning of the infall. There, repulsing forces such as 
the pressure barrier, the centrifugal force or the magnetic 
force are weaker than gravity. Therefore, Be starts with 
a negative value. When the gas gets closer, both the ki- 
netic and magnetic energies increase. Then Be becomes 
positive, which allows the gas to escape in the form of an 
outflow. It seems most natural that more powerful out- 
flows preferentially form along the rotational axis since 
there Be is most positive in all cases. Thus a bipolar mor- 
phology for outflows arises quite naturally in the present 
scenario. It is worth comparing this result with the zero- 
pressure limit ( Lery et al. 199S ) where Be is conserved on 



streamlines. In this case, the model does not produce ve- 
locities greater than the escape velocity anywhere. This is 
due to the fact that energy that is first gained by the mag- 
netic field at the expense of gravity and rotation during 
infall is subsequently returned as gravitational potential 
energy and kinetic energy in the outflow phase. 

The most rapidly outflowing gas is always near the 
symmetry axis because these streamlines pass closest to 
the star, deeper into the gravitational potential well. Also, 
the material on these streamlines is heated the most vigor- 
ously. We note that we do not include radiation pressure in 
our model. In the present context, the role of the radiation 
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Fig. 7. Gain in kinetic energy given by the ratio of 
the final to initial values in the flow, as a function of the 
initial angle of infall (left panel). Magnitudes of the kinetic 
energy at the turning point as a function of the initial 
angle of infall (right panel) (Same set of parameters as in 
Fig. |) 



produced by the source is to transport outward thermal 
energy. The radiation is eventually absorbed contributing 
to the gas heating. 

As the gas gets closer to the source, it rotates faster. 
The infalling plasma therefore has a larger electric current 
driven by the rotational motion. This increases the mag- 
netic energy, which is eventually converted into kinetic 
energy as the gas is redirected outwards. The magnetic 
field acts to collimate and accelerate the gas towards the 
polar regions. There the flow presents a strong poloidal 
velocity and a very low magnetic energy. 

In Fig. [?], we show the ratio of the final to initial mag- 
nitudes of the kinetic energy in the flow as a function of 
the initial angle of infall. The solution has been obtained 
for the same set of parameters as in Fig. |(| Since the gas 
falls from the equatorial regions and ends up in the axial 
zone, the ratio gives the gain in energy acquired during 
the circulation. The gain in kinetic energy varies drasti- 
cally. It goes from less than unity in the region far above 
the disk (0i n w 45°) to more than 10 for streamlines that 
start close to the equatorial region (9 in w 10°). In the lat- 
ter case, if the gas starts with an infalling velocity around 
5 km s , it will acquire an escaping velocity of at least 
50 km s _1 (|Bachillcr fc Tafalla 1999|). Again, the circula- 



tion is thus very efficient in accelerating the flow for the 
streamlines that start the closest to the equator, and then 
travel the closest to the central protostar, going deeper in 
the potential well. 

Therefore, in the axial region, the velocity will ap- 
pear to decrease as a function of the distance from 
the axis. Such properties has been observed for colli- 
mated molecular outflows (e.g. Guilloteau et al. 1997), 
and more recently, also for the central optical jet 



(Bacciotti et al. 200C). On the right panel of Fig. ^ the 
kinetic energy is represented at the turning point as a 
function of the initial angle of infall. The energy is maxi- 
mum for an intermediate range of angles, between 25° and 
35°. 

Summarizing, the material ejected in the outflow ac- 
quires a positive energy due to a combination of heating 
by the central protostar, pressure gradients developed by 
the infalling material, rotation and magnetic forces due 
to the non-zero Poynting flux. The study of the energetics 
along a streamline helps to describe how the local infalling 
flow is channeled into a bipolar outflow. 



5.3. The Electric Current and the Plasma j3 Parameter 

The electric current is a major indicator of the collimation 
of outflows. Indeed, it has been shown that the asymp- 
totic shape of steady axisymmetric magnetized outflows 
is either paraboloidal or cylindrical according to whether 
the electric current carried at infinity vanishes or not 
(See Lery, Heyvaerts, Appl & Norman 1999 and reference 
therein). In the present model, the axial electric current 
carried by the flow is given by 

I oc (r/r ) a+1/i , (19) 

with 1/4 > a > -1/2. If a < -1/4 (i.e., solutions with 
a gain in energy by radiation), the current is only large 
close to the protostar and will vanish at infinity. Thus the 
outflow should have a paraboloidal or conical shape. On 
the other hand, for the special case a — —1/4, i.e., solu- 
tions with radiative balance, the current is not a function 
of the distance to the central object but only a function 
of the angle. The current will have a finite value at in- 
finity and therefore the flow should focus cylindrically. A 
similar conclusion is obtained for a > —1/4 where the 
current increases with the distance. Consequently, for this 
range, the outflow should collimate cylindrically far from 
the source, at least in the axial regions. 

It is also interesting to study the plasma (3 parameter 
that corresponds to the gas to magnetic pressure ratio. 
For our model, (3 is given by 3 = y^P/u^. This quantity 
varies drastically with the angle in the circulation model 
as shown in Fig. [|. 3 starts from 10 in the infall region, 
reduces to 10 -2 at the turning point and finally increases 
up to several hundreds in the jet region. 

5.4. Down to the Thin Disk 

The last point that we shall discuss in this section con- 
cerns the boundary conditions at the interface with the 
thin disk, that in the model, is regarded as an open bound- 
ary. Physically, this corresponds to the entry of the gas 
into this region, while it is still falling onto the central 
object. The gas is not effectively falling onto the thin disk 
since self-gravity is not taken into account in the model. 
We have represented the shape of the poloidal component 
of the magnetic field above the disk in Fig. || for sev- 
eral solutions, that differ only by the initial conditions for 
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1 6.1. The Molecular Outflow Acceleration 

From an observational point of view, molecular outflows 
generally consist of a pair of oppositely-directed and 
poorly collimated lobes, which are symmetrically located 
about an embedded young stellar object, and detected 
in broad mm-wave emission lines, especially of carbon 
monoxide (CO) ( [Bachillcr fc Tafalla 1999) ). They are usu- 
ally moderately collimated (with initial opening angles 
9i « 20° — 90°), but highly collimated molecular outflows 
(8i < 10°) are also observed (e.g. Andre et al. 1990). The 
most luminous YSOs, which generally correspond to the 
most massive ones, have outflows with relatively wide- 



e 

Fig. 8. Variations of the plasma [3 parameter with 
the angle along an inner streamline starting from the 
cloud boundary for the circulation model. 6 m i n , 6 max and 
6t respectively correspond to the angles at the minimum, 
the maximum and the turning point of the solution (a = 
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Fig. 9. Shapes of poloidal sections of the mag- 
netic/flow surfaces above the disk for solutions of 
the infalling envelope model with different initial condi- 
tions for the the self-similar magnetic variable. The grey 
zone corresponds to the thin disk. In this solutions, the 
gas is flowing down as indicated by the arrow. (Same set 
of parameters as in Fig. |^) 

the the self-similar magnetic variable. We find that the 
angle of the trajectory of the infalling gas when it goes 
out of the boundary directly depends on the initial mag- 
netic field. For example, nearly perpendicular field lines 
can be obtained for a large magnitude of the field. When 
the magnetic field decreases, the lines bend and present a 
larger angle w.r.t. the perpendicular to the disk as shown 
in Fig. ||. In this case, the gas entering the disk has a non- 
vanishing radial velocity, and the magnetic field has an 
helical structure. 



6. Discussion 



opening angles (Churchwell 1997). Their outflow masses 
are estimated to be up to several hundred solar masses 
(e.g. 200M Q for Mon R2 (fTafalla et al. 1997j)), and in gen- 



eral they can be more than 5 times larger than the final 
mass of the central star. 

In the present model, most of the properties of molec- 
ular outflows can be understood in terms of the jet and 
the outflow having separate origins, although they may in- 
teract at their interfaces. Observationally, the spatial co- 
incidence of shocked molecular hydrogen bows with peaks 
in the CO outflow emission suggest that the jet and out- 
flow interact (e.g. Davis & Eisloffel 1995). In our model, 
the jet comes from the inner-most part of the central ob- 
ject, i.e., the inner accretion disk and the forming star, 
while the molecular outflow is produced by the circula- 
tion and is tied to the physics of the molecular cloud. We 
self-consistently model the infall, outflow, and jet regions 
using the same set of equations around the same central 
object. Even though the three different parts of the model 
are formally disconnected, they are related through strict 
boundary conditions at their interface surfaces. 

This has interesting consequences concerning the sub- 
sequent interaction of the flow with the coaxial jet. Since 
in our case, the difference in velocities between the jet and 
the molecular outflow material are reduced from the start, 
the shocks in the zone of acceleration due to their interac- 
tion should be less strong. The central fast jet has still the 
largest part of the total momentum per unit area, and the 
molecular outflow could undergo a prompt entrainment 
from the head of the jet. The effects of this entrainment, 
however, are reduced because of the motion induced in the 
CO outflow by the circulation around the protostar. 

But the most interesting feature of the circulation 
model is probably that it could produce solutions where 
the mass of the molecular outflow is larger than the final 
mass of the forming star. This would particularly true if 
self-gravity were included in the model. One might then 
understand how bipolar outflows from massive protostars 
are observed to transport masses largely exceeding those 
of the associated stars. For high luminosity objects, so- 
lutions can be obtained that show well-collimated jets, 
even though the CO flow may appear poorly collimated 
QBachiller fc Tafalla 19991 ) . 

Other important features of the model are that the 
distribution of the flow mass with velocity approximately 
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follows a power law, and the flow collimation increases 
with velocity. The model also shows solutions in which 
the flow velocity increases with distance along a line par- 
allel to the outflow axis. Also, the synthetic CO spec- 
tra corresponding to these solutions nicely reproduce 
most of the observational features of molecular outflows 



(Fiege & Henriksen 1996, Lery et al. 1999). The acceler- 
ation is relatively gentle in this type of model, so that 
molecules are advected with the flow with relatively little 
destruction. This is in contrast to the more severe accel- 
eration that occurs when gas is entr ained by a fast jet , 
which tends to destro y the molecules ( Tafalla et al. 1997 , 
Downes fc Ray 1999 ). In our model, several tens of solar 
masses are moving at supersonic velocities while remain- 
ing neutral and in molecular form. The calculations pro- 
duce solutions where the outflow can have large opening 
angles, and where the most massive protostars produce the 
fastest and the less collimated outflows ( Lery et al. 1999| ). 

In conclusion, we suggest that CO outflows are domi- 
nated by the global circulation of material around the pro- 
tostar, except for in a thin layer surrounding the jet, where 
the dynamics is governed by entrainment. We stress that, 
in the present model, the two flows (atomic and molecu- 
lar) are not strongly linked dynamically, and, hence, there 
is no need to transfer large momentum from the jet to the 
molecular outflow through the entrainment processes. 



6.2. Molecular Cavities 

At CO mm- wavelengths, e.g. L1551 

( [Moriarty-Schieven fe Snell 1988|, |Hartigan et al. 2000D , 



NGC 2071 (Moriarty-Schieven et al. 1989), or Mon 



R2 (Tafalla et al. 1997 ) and also in the infrared band 
flClose et al. 1997a|) , cavities inside molecular outflows 
are observed. The edges of such cavities, located between 
the jet and the molecular outflow, are delineated by 
H2 emission. Indeed, close to the source, "throats" 
are o bserved 100 to 200 AU above the accretion 
disk ( Close et al. 1997b ). The association of outflows 
with cavities is not limited to low-luminosity sources 
( ICampbell et al. 1986j [Lafalla et al. 1997|). For exa mple, 



the blue lobe of AFGL 490 ( |Campbcll et al. 1986[ ) coin- 
cides with a conical nebula of scattered light, suggesting 
that the outflow interacts with a shell of dense gas 
surrounding a partly evacuated channel. 

The cavities can have large opening angle in the range 
90°-125° (e.g. Velusamy fe Langer 1998), and their na- 
ture is unclear ( Chandler et al. 1996| , pomez et al. 1997 



Gucth et al. 1998). These observations are in good agree- 



ment with the opening angles in our model and in wide- 
angle wind models, which produce wider cavities than 
jet models. It has also been proposed that the cavities 
could be produced by the wandering of a precessing jet. 
However, it appears that the dynamical time-scale of the 
jet is too short for this explanation to work in some ob- 
served sources (IChandler et al. 1996|) . 



In the present model, the molecular cavities may be 
identified with the substantial decrease in density in the 
intermediate region between the jet and the molecular out- 
flow (See Fig.|| in the axial region). Note that solutions 
with wider opening angles than the one presented here can 
easily be obtained. Moreover, the central jet region may 
be comprised mainly of atomic gas, which is not observ- 
able in molecular lines. This is likely because the gas will 
probably be dissociated and ionized near the protostar. 
As the gas escapes and cools, atomic gas will form rapidly. 
However, the timescale for the production of molecules is 
much longer. Thus, the molecular outflo ws appear as a 
hollow conical structure ( Lery et al. 1999 ). The model in- 
dicates a wide range of opening angles from 40° to 140°, 
depending on the input parameters, in agreement with 
the observed values. However, this angle may vary with 
time but the timescale for the widening is expected to be 
much longer than the typical dynamical timescale. Thus 
our model may represent a quasi-steady flow pattern at 
some instant in a slowly changing outflow. Anyhow, the 
cavity is a result of the circulation itself and we do not 
need to assume that the jet is precessing. 



6.3. Asymmetric Outflows 



Another important observational issue is the large num- 
ber of bipolar optical outflows showing properties that dif- 
fer substantially in the blue-shifted and red-shifted parts 
(Hirth et al. 1994, Bacciotti et al. 1999). In the axisym- 



metric quadrupolar circulation model, solutions are not 
limited to an hemisphere, and consequently can differ 



from side to sid e of the disk ( [Fiege fc Henriksen 1996 



Lery et al. 1999|) . . We present an example of such a solu- 
tion in Fig. |lO|. In this figure, the streamlines are rep- 
resented for all the three different parts of the model, 
i.e., the circulation model, the jet model and the infall 
envelope. The solutions have the same set of parameters 
on both sides of the equator, i.e., a = — 0.12,9 = 0.35, 
r Q ss 100AU, M ss 1M@, but have different opening 
angles, that we prescribe as input parameters (6 max — 
27r/3 > ir/2 for the upper part, and max = ir/3 < ir/2 
for the lower part). These solutions may be generated by 
different conditions on both sides of the disk, for exam- 
ple in a region of a cloud that had a strong pre-existing 
density gradient, or at the edge of a molecular cloud. 
The latter case is similar to the standard interpretation 
of the asymmetric objects as an outflow breaking out of 
the cloud on one side. Note, however, that contrary to 
the optical jet case, for the molecular outflow the weak 
lobe will be pointing out of the cloud as observed for 
many unipolar CO outflows known, such as NGC2024- 
FIR5 (|Richer et al. 1992|), or the al most-unipolar HH46- 
47 system ( Chernin fc Masson 1991 ). Our model also pre- 
dicts that the dominant outflow lobe will be on the side 
pointing into the cloud because that is the side where ma- 
terial can be gathered. Thus, the various asymmetries ob- 
served on opposite sides of many bipolar outflows, could 
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Fig. 10. Streamlines of solutions corresponding to 
asymmetric outflows with quadrupolar geometry. The 
parameters are similar to the previous solutions (a = 
— 0.12, Oo = 0.5) but the maximum opening angles of the 
solutions are not limited to 90°. 



find a natural explanation in the quadrupolar geometry. In 
any event, detailed comparisons of solutions for a specific 
asymmetric source could be a good test of our hypothesis. 
It is also noteworthy that purely unipolar models had al- 
ready been found in the context of the circulation model 
when treated without Poynting flux (Ficge 1996). 



6.4. Collimation 

Several mechanisms, from purely hydrodynamic 
(Mellema fc Frank 1997) to magnetic ones, have 



been proposed to account for the collimation of jets. 
Magnetic processes seem very promising and they 
are generally accepted as the primary option. In this 
case, the outflow is collimated by the pinch effect 
of the toroidal component of the magnetic field that 
arises as a consequence of the rotation of the disk 
(Lery et al. 1999, Breitmoser & Camenzind 2000). In 



cylindrical jet asymptotics, the outflowing plasma reaches 
an equilibrium wherein it is confined by magnetic forces 
or gas pressure gradients, while it is supported by cen- 
trifugal forces or gas pressure gradients. The idea that the 
external medium has, in fact, a non- negligible effect has 
been raised recently (Okamoto 1999, Lery fc Frank 2000| ). 
Indeed, magnetized winds may not cylindrically collimate 
without an external help, such as the channeling effects of 
thick accretion discs and/or confinement resulting from 
the ambient medium. Moreover, in spite of the success 
of magnetic models, recent numerical studies have shown 



that pure hydrodynamic collimation could be effective at 
producing jets ( [Icke et al. 1992[ , pVlcllcma fc Frank 1997| ), 
thanks to the "shock focusing" mechanism due to the 
interaction of the outflow with the external environment. 
In the present model, the collimation of the axial jet 
arises to a large extent from the hoop stress of the 
toroidal field component but it is also reinforced by the 
density gradients in the surrounding molecular outflow. 
Thus, both hydrodynamic and magnetic effects play 
a role. Note that the present model predicts that the 
collimation does tend to increase with the velocity of 
t he gas regardless of the mass of the central object 
flFicgc fc Hcnrikscn 199l| |Lery et al. 1999|) . 

The presence of low density polar regions may also 
give a preferential direction and help the alignment of jets 
coming from a binary system, as in the case of L1551. 
Here, HST imaging and spectroscopy from the ground 
( Fridhmd fc Liseau 1998 , Reipurth et al. 200C ) indicates 
the presence of two individual jets almost aligned (an- 
gle between the jets less than 20°) and lying inside a CO 
cavity. 



6.5. Possible origin of the quadrupolar topology 

Most of current outflow models are based on a dipolar ge- 
ometry where the magnetic field lines are either advected 
via accretion towards the central object in the accretion 
disk or are created in situ. Aburihan et al. (2001) proposed 
a simple, qualitative mechanism to generate a quadrupo- 
lar field topology, which we briefly summarize here. The 
basic idea is that advection deforms the field into a pre- 
dominantly quadrupolar topology as the magnetized gas 
collapses. During the earliest stages of formation, the mag- 
netic field near the protostar would be mainly dipolar 
as the gas starts to collapse. Since the magnetic field is 
frozen into the gas, the infalling material advects the mag- 
netic field toward the accreting core until it is strongly 
pinched toward the protostar in the equatorial region. At 
this point, the field is mainly radial. Poloidal pressure gra- 
dients, with both radial and angular components, develop 
as the protostar forms. Eventually these pressure gradients 
become strong enough to deflect some of the gas toward 



the polar regions (Lery et al. 1999). The magnetic field 
near the protostar is advected with this outflowing gas, 
and eventually forms magnetic arches near the polar axis. 
A predominantly quadrupolar field topology remains as 
these magnetic arches are advected to large radii. Other 
mechanisms might exist to convert an initially dipolar field 
into a quadrupolar geometry, but these might rely on re- 
connection or dynamos. 

7. Conclusion 

We have proposed a new global MHD model for flows 
around YSOs. The global model combines a jet model, a 
circulation model and an infalling envelope model. Instead 
of the usual mechanisms invoked for the origin of molecu- 
lar outflows, the outflow is powered by the infalling matter 
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through a heated quadrupolar circulation pattern around 
the central object. The molecular outflow may still be af- 
fected by cntrainment from the fast jet, but this would 
be limited to the polar regions and it would not be the 
dominant factor behind its acceleration. 

The solutions are developed within the context of r- 
self-similarity. We use the set of steady axisymmetric 
MHD equations for all the three parts of the model. The 
only physical scales that enter into our calculation are 
the gravitational constant G, the fixed central mass M, 
and a fiducial radius r a . The parameters of the model are 
the indexes a and a/ of the self-similar system, and m i n 
and 9 max , the limiting angles between the three regions 
of the model. Magnetic field and streamlines are required 
to be quadrupolar in the poloidal plane for the circulation 
model. 

We have explored the physical behaviour of the model 
by studying the flow through its various energy compo- 
nents along streamlines. We have reported here solutions 
for the three different parts of this self-similar model. An 
example of solution is presented for a one solar mass ob- 
ject, even though this class of model can apply to both low 
and high mass protostars. The solutions show dynamically 
significant density gradients in the axial region, precisely 
where the radial velocity and collimation are the largest. 
Rotation is maximum at the turning point but it is still 
dynamically important close to the disk. 

Our model can reproduce mass loss rates that arc 
in agreement with values derived from observations. 
For the specific case that is presented here, we obtain 
Mj e t/Macc ~ 0.1, as observed in several systems, even 
if our model does not deal with the precise jet launching 
mechanism. Other values, typically ranging from ~ 0.04 
to ~ 0.4 are possible, which may impose additional obser- 
vational constraints on the parameter space of our model. 

From an observational point of view, we clarify the 
nature of the molecular outflow acceleration and its rela- 
tion with the fast jet, by providing a global picture of the 
jet/outflow system, which does not primarily rely on cn- 
trainment (prompt or turbulent). In this scenario, the ob- 
served molecular cavities result from the circulation pat- 
tern itself, and there is no need for jet wandering. Our 
model helps to explain the presence of high mass molecu- 
lar outflows around massive objects, being the molecular 
material simply circulating around the central object from 
the nebula to the bipolar outflows. The various asymme- 
tries observed on opposite sides of many bipolar outflows 
could find a natural explanation in the quadrupolar geom- 
etry; our model predicts that the dominant outflow lobe 
will be on the side pointing into the cloud. The present 
study suggest that radiative heating and the Poynting 
flux may ultimately be the main energy sources driving 
molecular outflow. Finally, although the details of the jet 
mechanism may be peculiar to individual objects we be- 
lieve the infall-outflow circulation to arise naturally given 
accretion, and thus could also be present in AGN's and 
around suitably placed objects such as SS433. 
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